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We present a solution to the problem of broadband decoupling of a coupled homonuclear two-spin 
system. We describe a pulse sequence that creates an effective field perpendicular to the coupling 
interaction with a magnitude propotional to the chemical shift of the spins over a broad range of 
chemical shifts. When the chemical shift difference as imprinted on the perpendicular field is greater 
than the coupling between the spins, we get effective decoupling. The pulse sequence may be useful 
in various NMR applications. 
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Coupled spin systems are ubiquitous in nuclear mag- 
netic resonance (NMR) [HE]. In NMR, the coupling be- 
tween spins causes multiplcts in the spectrum of the Free 
Induction Decay (FID) signal [2J, which complicates the 
interpretation of the data. For example to determine the 
structure of a large protein, it is desirable that these mul- 
tiplets be collapsed into singlets, thereby increasing the 
resolution and the intensity of the peaks. 

If the coupled spins belong to different spin species 
(heteronuclear spins), their Larmor frequencies are well 
separated. The underlying concept in many commonly 
used decoupling techniques [3HZ] is inverting one spin 
rapidly while observing the other spin. Because the Lar- 
mor frequencies are very different, the radio frequency 
(rf ) radiation can invert one spin species without affect- 
ing the other 

Larmor frequencies for Homonuclear spins are only 
separated by the chemical shifts, which are caused by 
different chemical environments around each spin that 
alters the strength of the magnetic field at their loca- 
tions. Even at high magnetic fields used in the state of 
the art NMR experiments, the smallest chemical shift 
difference can be of the same order as the strength of the 
couplings between the spins. Therefore, it is very hard to 
selectively invert one spin without disturbing its neigh- 
bors, especially when the precise chemical shifts are un- 
known and are the parameters to be determined. There 
is some success in the use of band-selective shaped pulses 
for decoupling of a coupled IS homonuclear spin system, 
where spin / belongs to a band of frequencies which is 
well separated from spin S that lies in a distinct band 
(see [8], [9], [10]). This may be the case in many practi- 
cal applications. However, when there is no information 
about the location of the chemical shifts, it is not clear 
how to do the selective pulses. Here, assuming almost no 
knowledge about the chemical shifts, we present a new 
pulse sequence that decouples homonuclear spins with- 
out the need for selective inversion over a broad range of 
frequencies. 

The Hamiltonian H a of a system of two spins I and 
S coupled together in the rotating frame at the Larmor 



frequency of the spin species is 

H = ljjI z + u> s S z + 2irJI z S z (1) 

where ujj and u>s are the chemical shifts of the spins / 
and S. 1/3 and Sp, f3 = x,y,z are the spin operators. 
J is the strength of the scalar coupling. This coupling, 
also named indirect interaction due to shared electrons 
in the bonding of diffrcnt atoms, has the general form 
2ttJ(I • S) = 2nJ(I x S x + I y S y + I Z S Z ). However, under 
the condition 

2ttJ « \uj - u>s\ (2) 

it is truncated to 2irJI z S z . This is called high-field, or 
weak-coupling approximation For this paper we as- 
sume that we have the Ising coupling as in Eq. ([I]) and 
present a pulse sequence that decouples this coupling as- 
suming no knowledge of the chemical shifts, except the 
relation Q and that the chemical shifts fall in a certain 
range. We then show that even in the presence of general 
isotropic interaction 2nJ(l ■ S), we can create an effec- 
tive zz interaction to which our pulse sequence can be 
applied. 

We will attempt to create an effective Hamiltonian 

H^t) = -(ujjIy+ujsSy) + 2nJI z S z . (3) 

This Hamiltonian represents an effective y field which is 
different on two spins propotional to the chemical shift 
difference and perpendicular to the zz coupling Hamil- 
tonian. If we transform to the interaction frame of this 
field, we average the coupling. To synthesize the Hamil- 
tonian Hi, we evolve rf Hamiltonian A(I X + S x ) for first 
2 At units of time and — A{I X + S x ) for next 2 At units 
of time, where At is appropriately chosen. We make 
the chemical shift Hamiltonian evolve as wil z + wsS z , 
—u)il z - u>sS z , - usS z and UtI z + ojsS z , for At 

units of time each. We now observe how the rf Hamil- 
tonian evolves in the interaction frame of the chemical 
shift Hamiltonian. This can be seen from the Fig. [TJi. 
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Assuming A At is a small flip angle, the x component of 
the rf Hamiltonian cancels while the y component adds 
at the end of the cycle. If a = uiiAt is a small tip angle 
then the y component evolves as ~ Aujjt which over a pe- 
riod At integrates to ^^uJiAt = |cjjAt. This produces 
an effective Hamiltonian + tu B S v ). In the inter- 

action frame of the chemical shift, the coupling terms 
by themselves simply do not evolve. This then creates a 
field perpendicular to the couplings carrying the chemical 
shift difference which helps to decouple the spins. Chem- 
ical shifts can be inverted by simply applying a hard tt 
pulse. This doesnot affect the coupling Hamiltonian. Fig. 
[TJj shows the basic pulse sequence with n pulses with x 
phase that inverts the chemical shift Hamiltonian and 
produces the cycle as described above. Similar cycles of 
the rf Hamiltonian have been recently proposed to de- 
sign small tip angle-excitation pulses that are robust to 
rf inhomogeneity [T2] . 
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FIG. 1: la shows the evolution of the rf Hamiltonian (only 
I component is shown) in the toggling frame of the chemical 
shift Hamiltonian that goes through a cycle (ojiI z + u>sS z ), 
-~{ujiI z + uj s S z ), —(uil z + usS z ) and (uih + ojsS z ), for an 
appropriate period of At units each. The angle a = uiiAt. 
lb shows the basic pulse sequence that makes the chemical 
shift and rf Hamiltonian excecute the cycle in Fig. la. The 
basic building block is repeated before the precession of the 
spins and the resulting chemical shift is sampled and the unit 
is repeated. 

Assuming chemical shifts are dispersed over the range 
[-B, B], i.e, \vj\, \ug\ < B, we now compute the effective 
Hamiltonian resulting from the above pulse sequence. We 
evolve the system under the following four Hamiltonians, 
each in a small time interval At. 

iJi =mh + u s S z + 2ttJI z S z + A(I X + S x ) 
H 2 =- uih - u s S, + 2ttJI z S z + A(I X + S x ) 
H 3 = - ujl z - uj s S z + 2irJI z S z - A(I X + S x ) 
Hi =uj t I z + lo s S z + 2nJI z S z - A(I X + S x ) 



(4) 



The unitary operator generated by the four period 
pulse sequence above is 



U 



-iH 4 At -iH 3 At -iH 2 At -zHi.At 



(5) 



We wish to find the effective Hamiltonian 



U = e - l4 ( At ) ff eft 



(6) 



Using the Baker Campbell Hausdorff formula and 
keeping terms up to the third order in At, where 
VB 2 + A 2 At < 1, we have 

U = CTq>{-4i(At)2TrJI z S z - 2i{At) 2 A( Wl I y + w s S v ) 

- 4i{At) 2 A2TrJ(I y S z + I z S y ) + 2i(At) 3 A 2 { Wl I z + w s S z ) 

- l6i{At) 3 A 2 2TTj(I y S y - I Z S Z ) + 0(A£) 4 ) 



(7) 



Compare the above equations, we have 
1 



H. 



eff : 



2ttJI z S z 



-MUlly 



e2TTj(i y s. 
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IzSy 



+ U S S y ) 

< 2 (^I z +oj s S z ) (8) 



2 2TTj{I y Sy ~ I Z S Z ) 



where as before 9 = AtA. Ideally, we only want the 
first two terms: the J coupling stays in the z direction, 
while the main chemical shifts being scaled down by \Q 
effectively point in the y direction. There is a residual 
chemical shift still remaining along the z direction and 
therefore we get an effective field in the direction y' which 
makes an angle 7 to the y axis. We resolve the coupling 
parallel and perpendicular to this effective field. The per- 
pendicular part is truncated and the parallel component 
reduces the magnitude of residual coupling from ^82ttJ, 
to ^82irJ as described below. We write the full Zeeman 
Hamiltonian as ^(wjlyi +wsS y i), where we neglect small 
change in 9 arising from adding 9 2 in quadrature. Then 
ly = ly' cos 7 + I z > sin 7, I z — I z i cos 7 — Iyi sin 7 where 

kg 2 

tan 7 = ^-a = 9. Replacing y and z in terms of y' and 
2 H 

z', the coupling term in the y direction becomes 

= J sin 2 7 — 92tt J(sin 7 cos 7 + cos 7 sin 7) 

+ -9 2 2n J(cos 2 7 - sin 2 7) 
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= 2ttJ(1 
+ ^6» 2 J(2cos 2 7 



2ttJ 



cos 2 7) — 292ttJ tan 7 cos 7 
-1) 
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J J(l-2 + 4/3) 
1 + 9 2 
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So the most important terms in the effective Hamiltonian 
are 



H't. 



+ -(iU I Iy> +UJ S Sy>) + —2-KjIy,Syl 



(9) 
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That means, the chemical shifts are scaled down as 9, 
while the coupling strength is scaled down as 9 2 . This, 
as shown in the simulation below is a big improvement 
in enhancing the spectrum. The smaller the value of 9 
, the more effective the decoupling is. Scaling down the 
chemical shifts too much is not preferable as it degrades 
the resolution. Note that the chemical shift information 
of the spins is still retained and since the tilt 9 is well 
defined (assuming no rf-inhomogeneity) we can make the 
effective Zeeman interaction point along z axis by appro- 
priate rf rotations. 

Figures 2 and 3 show the results obtained by simulating 
the decoupling sequence starting with an initial density 
operator p(0) = I x + S x . We simulate the evolution of 
• S W = (-^x(^) + S x (t)). The simulation parameters have 
been chosen to reflect true experimental conditions. In 
absence of decoupling the original signal has an envelope 
of cos(nJt), that of the decoupled signal decays much 
more slowly: approximately cos[^(AAt) 2 nJt}, which is 
our analytical calculation and is 80 times as long as the 
original one. The spectrum (Fourier transform of s(t)) 
of the decoupled system shows sharp peaks with much 
higher intensities. 

Notice that, in the spectrum of the decoupled system, 
we choose to truncate the signal after time 20 j. This 
is a very good assumption because in practice, the FID 
is only collected upto the transverse relaxation time T 2 , 
and 20 j >> for the coupling values and relaxation 
rates of interest. 

In summary, we have seen that as the Larmor frequen- 
cies are scaled down by a factor of ^(At)A, the coupling 
strength J is scaled down by a factor of ^(At) 2 A 2 . For 
protons, chemical shift is 10 ppm, so for a magnetic field 
B Q = 500MHz, the Larmor frequencies are spread over 
a range B — 5kHz. The coupling strength is of order 
J 50Hz. In NMR, the static magnetic field can be 
made homogenous to 1 part per 10 9 , so with a magnetic 
field of 500MHz, we can resolve peaks that are up to 
O.bHz apart. If ^((At)A) 2 is then the signal envelop 
falls down in order of a second which is enough to resolve 
frequencies few Hz apart. This is satisfied by the choice 
of (At) A — |. In the simulation, we choose (At) A = |, 
which is very close to the limit here. We assume B = 1 
kHz. We choose A = B. Smaller values will not satisfy 
VA 2 + B 2 At < 1 in Eq. Larger value of A is not 
preferable as it gives smaller value of At , for a fixed 9, 
which means more building blocks and additional ir pulse 
overhead in the blocks. If we do the analysis of effective 
Hamiltonian in the toggling frame of the Chemical-shifts, 
then we can choose AAt as an order parameter and still 
have BAt > 1. This however gives a frequency depen- 
dent scaling of the chemical shift. 

In general, the coupling Hamiltonian between 
Homonuclear spins is the isotropic Hamiltonian IS. If we 
let this Hamiltonian evolve freely under chemical shifts, 
the interactions can be averaged to Ising interactions. We 
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FIG. 2: (a) Signal (time in units of j), (b) Spectrum without 
decoupling taken over time 4. U)i = 27rl20,a;s = 27rl00, J = 
1 (rescale time to normalize J = 1). 



can combine two such free precession periods of length t\ 
and T2 such that t± — t 2 = At by placing a ir pulse in 
between as shown in Fig. |4^. The x-pulse at the end of 
the period with flip angle 9 produces a net Hamiltonian 
~ {ujI.+u.S, + ^2irJI z S z + A(I X + S x ) + m (I X S X + 
lySy) + K2(I x S y — I y S x ))At which is the basic building 
block of our pulse sequence. The couplings are scaled by 
a factor Tl ^J 2 , Ki and k 2 depend on ujj — uj s and time T\ 
and T2 as shown in[4j;. By choice of different t[s and r^s, 
over various building blocks, we can make k% ~ over a 
broad range of Wj — u>s- Similarly we can make k-i ~ 
(or invert the sign of «2 with a y phase pulse to reduce 
scaling of the zz couplings). 

The decoupling sequence scales down the chemical 
shifts by a factor 9. This obviously degrades the resolu- 
tion but because we can decouple effectively we can make 
spins precess for longer periods and restore resolution. 
The most ideal setting to do this is to transfer the spin 
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FIG. 3: Results using decoupling sequence wi = 27rl20, ujs = 
2tt100, J = 1, A = 10 3 , At = 2 x 10~ 4 . (a) Signal that de- 
cays much more slowly with the envelope of cos(|(AAi) 2 7r Jt), 
(time in units of j) (b) Spectrum obtained by truncating the 
signal at time . n pulses are assumed to be delta pulses. 



where the true spectrum comprises of spikes at N(f) 
is the noise spectrum. We wish to deconvolve the ob- 
served spectrum Z(f) to recover X(fi). To do this, first 
calibrate %(/) from an isolated peak. x(f) reflects the 
inhomogeneity of the sample. We can then deconvolve 
using least square estimator. 

In conclusion, we have introduced a new pulse se- 
quence to decouple homonuclear spins with Ising inter- 
actions and proposed how the solution can be extended 
to Isotropic couplings. 
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FIG. 4: Panel a shows how an isotropic coupling Hamiltonian 
may be averaged to an Ising coupling by free precession pe- 
riods, where the chemical shift difference between the spins 
evolves the planar part xx + yy of the isotropic iteraction. 
Panel c shows the evolution of the planar Hamiltonian xx+yy 
as result of the building block in Panel a. Panel b shows the 
application of the decoupling sequence in the indirect evolu- 
tion period, where spins precess for constant time. The center 
7T pulse refocusses the chemical shift so there is chemical shift 
evolution for time di — d,2- 



coherence after it acquires a phase under chemical shift 
evolution on to a coupled spin for readout, as in the in- 
direct evolution periods in multi-dimensional NMR. The 
basic setup is as shown in Fig. where spins precess for 
constant time d\ + d^. Moving the center tt pulse causes 
chemical shift evolution as function of difference d\ — di. 
Since it is a constant time evolution , the relaxation loss 
is constant and hence does not effect resolution. Effective 
decoupling ensures that d\ + di can be made large. 

Since the chemical shifts are scaled by a factor 9 = 
AAt which depends on the strength of the applied rf 
field, in the presence of rf inhomogeneity, the observed 
spectrum is blurred by the point spread function x(/). 
The observed spectrum is therefore 

Z(f) = J2x(fiMf-fi) + N(f) (10) 
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